Low-energy properties and magnetization plateaus in a 2-leg mixed spin ladder 
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Using the density matrix renormalization group technique we investigate the low-energy properties 
and the magnetization plateau behavior in a 2-leg mixed spin ladder consisting of a spin- 1/2 chain 
coupled with a spin-1 chain. The calculated results show that the system is in the same universality 
class as the spin-3/2 chain when the interchain coupling is strongly ferromagnetic, but the similarity 
between the two systems is less clear under other coupling conditions. We have identified two types 
of magnetization plateau phases. The calculation of the magnetization distribution on the spin- 
1/2 and the spin-1 chains on the ladder shows that one plateau phase is related to the partially 
magnetized valence-bond-solid state, and the other plateau state contains strongly coupled S=l and 
s=l/2 spins on the rung. 



PACS Numbers: 75.10.Jm, 75.40.Mg 
I. INTRODUCTION 

The field of low dimensional quantum spin systems 
has been a focus in condensed matter physics for about 
two decades since Haldane pointed out the difference 
between the integer spin Heisenberg antiferromagnetic 
(HAF) chains and the half-integer spin chains n By map- 
ping the HAF chain to a non-linear a modelJH Haldane 
conjectured that there is a finite gap between the ground 
state and the low-lying excitation states for Heisenberg 
chains with integer spin, while it is gapless for half-integer 
spin chains. This conjecture has been verified by later 
analytic, numerical and experimental studies. _ 

In recent years, the physics of spin laddersB consist- 
ing of different number of coupled HAF spin- 1/2 chains 
has attracted much attention. The spin ladders of even 
number coupled chains have a finite gap in low-energy 
spectrum and those of odd number coupled chains are 
gapless. The comparison of the quantum spin chains 
with the spin ladders suggests that there may be some 
relations between the two systems. For the 2-leg lad- 
ders, White showed numerically that its ground state 
can transform continuously to that of the S=l spin chain 
no matter the interchain coupling of the ladder is ferro- 
magnetic or antiferromagneticjj therefore they are in the 
same phase. While by considering a more general frus- 
trated 2-leg ladder case, Wang identified two different 
phases, the Haldane phase and the singlet phase, cor- 
responding the ferromagnetic awt,the antiferromagnetic 
interchain coupling respectivelyu'tl 

The effect of the magnetic field on the spin chains and 
the spin ladders has been a topic of great interest. In 
the presence of a magnetic field, a spin- 1/2 chain will not 
open a gap until the applied magnetic field is greater than 
the saturate field. While for a spin-1 chain, because of the 
existence of the Haldane gap, there is a critical field where 
the gap is closed and the system remains gapless from the 



critical field to the saturate field. Below the critical field, 
there is a plateau at magnetization per site m=0. For 
those chains with higher spins, bond-alternating or other 
additional interactions, and the spin ladders, the mag- 
netization procesSr«.ay also display plateaus at non-zero 
m.Ena It is shownliJ that for a spin chain with periodic 
structures, in a uniform magnetic field the magnetization 
plateaus may exist at the magnetization per site m when 
the condition 
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is satisfied. Here n is the period of the ground state and 
S is the site spin. For ladder systems, the condition (|l|) 
needs some slight modification.^ The simplest case where 
a non-zero m magnetization plateau may exist is the tra- 
ditional translationally invariant S=3/2 spin chain, with 
n=l, and m=l/2 satisfying Eq. (|l|). But detailed anak 
yses have shown that there is no plateau in this caseE2l 
except when single-ion anisotropy D is present and D > 
0.93JD 

Considerable attention has been given to mixed spin 
chains in recent years. Materials with-apixed (Si,S 2 ) 
chain structures have been synthesized. 1111113 Both types 
of excitations, gaped and gapless, in the low-energy 




lternating spin 
Magnetization 
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spectrum of the fcrrimagnetic quanj. 
(1,1/2) chain have been obtainec 
plateau behavior in such systemsOi 
ternating spin chains also has been studied.EHl 

In this paper, we present the results of the low-energy 
properties and the magnetization plateau behavior of a 
mixed spin ladder consisting of a spin- 1/2 chain coupled 
with a spin-1 chain 
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using the density matrix renormalization group (DMRG) 
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method.0 Here S and s denote the spin operators on the 
spin-1 and spin-1/2 chains. ./V is the number of transla- 
tional invariant cell (S,s)=(l,l/2). The number of total 
lattice sites is L = 2N. We consider the case of antifer- 
romagnetic intrachain coupling J\, Ji > in this work. 
To simplify the discussion, we fix J\ + Ji—1 and make 
Ji change its value from to 2. The interchain coupling 
J 1 can be ferromagnetic or antiferromagnetic. It is of 
considerable interest to examine the effect of the nature 
of the interchain coupling on the properties of the mixed 
spin ladder, especially the possibility of introducing new 
phases as in the case of the spin-1/2 isotropic ladder .atl 
The mixed spin ladder is illustrated in Fig. [I]. From 
intuitive point of view, when the interchain coupling is 
strongly ferromagnetic this mixed spin ladder is expected 
to be equivalent to the spin-3/2 Heisenberg chain, which 
is gapless in the low-energy spectrum with a spin wave 
velocity v =3.87 and belongs to the same universality 
Haldane class as the S=l/2 Heisenberg chahi.til Our cal- 
culated low-energy properties provide strong support for 
this view. When the interchain coupling is antiferromag- 
netic, the calculated low-energy spectrum does not show 
clear evidence for the equivalence. 



plateau behavior and the phase diagram are discussed in 
Sec. [H. Finally, in Sec. [V, we present the conclusion 
of the calculated results and some discussion on the 2-leg 
mixed spin ladders with general mixed spin cell (S, s). 



II. LOW ENERGY PROPERTIES 

The ground state and low-energy excitation states of 
Hamiltonian (||) are calculated using DMRG by keeping 
400 states. We use the open boundary condition (OBC) 
and calculate up to the cell length N=60. The largest 
truncation error is of the order of 10 -7 . 

We calculated the energies and correlation functions 
of the lowest state in Sj: ot =0, 1, 2, 3 sectors for different 
Ji, J2 and J'. For all the parameters we calculated, 
the low-energy spectrum is gapless at the thermodynamic 
limit and except for those very small J', it is similar for 
different J%, and J2. In the following we present the 
calculated results of the excitation energy of low-lying 
states, the spin excitation in the system, and the spin- 
spin correlation functions for the 2-leg mixed spin ladder 
with Ji = J2 = 1. 
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FIG. 1. The illustration of the 2-leg mixed spin ladder. 
Large (small) circles denote S=l (s=l/2) spin. Intrachain 
coupling Ji, J 2 > 0. 

The magnetic plateau study shows that the J' >0 and 
J' <0 phases behave differently. In the presence of a 
uniform magnetic field H, an additional term 



H' = -hJ2(S! + sl) 



(3) 



needs to be added to Hamiltonian (g). When Ji and 
J' equal zero, there is a m=l/2 magnetization plateau 
from magnetic field H=0 to the critical field of the 
spin-1 chain, which is related to the partially magne- 
tized Kaleneei-bond-solid (VBS) state of the S=3/2 spin 
chain .II2E3c-3 In our calculation, besides this magnetiza- 
tion plateau phase, we also have identified another mag- 
netization plateau phase when the interchain coupling is 
strongly antiferromagnetic. We have calculated the mag- 
netization distribution on the two chains at m = 1/2. In 
the first plateau phase, the magnetization is mainly from 
the spin-1/2 chain as expected. But in the second plateau 
phase, the magnetization comes mainly from the spin-1 
chain, while the spins are antiparallel to the applied field 
in the spin-1/2 chain. 

We will show the calculated low-energy properties of 
the 2-leg mixed spin ladder in Sec. 0. The magnetization 
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FIG. 2. Excitation energy times the number of cells 
N(E n — Eq) for n=l (filled circle), n=2 (open circle) and 
n=3 (open square) are plotted versus 1/lnTV for Ji = 1.0 and 
(a) J'=-4, (b) J' =4. 

The excitation energy of several lowest states times cell 
length N(E n - E Q ) as a function of 1/lniV for J' = -4 
and 4 are shown in Fig. |^(a) and (b) respectively. Here 
the energy of the lowest state for Sf ot = k sector is de- 
noted as Ef.- For J' — —4, therfaswlts are nearly identical 
to that of the spin-3/2 chain. E3'c!3 Because we use OBC 
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in the calculation, there is a s=l/2 free end spi 
can be understood by valence bond solid pictureE 
The S* ot =0 and S* ot =l state will be degenerate in the 
thermodynamic limit. The energy spacing of the two 
quasi-degenerate ground states scales as l/(iVln_AQ, ex- 
actly same as that for a pure spin-3/2 chain.E3cJ For 
J'=4, there is no end spin. Here the S^ ot =0 state is the 
only ground state in the thermodynamic limit, and the 
first excitation state at finite chain length is a true exci- 
tation state. The energy spacing between this state and 
the ground state no longer scales as l/(iVln7V). 

The difference of the expectation value of the z com- 
ponent of each spin site in the excitation states to the 
ground state is shown in Fig. ||. For J'=4, due to the 
reflection symmetry of the system, the expectation value 
(Sf) for Sf ot = state is zero for all i, we show (Sf) di- 
rectly. For J' = —4, considering the existence of the end 
state, we show (Sf) n — {Sf)i. These differences exhibit 
spin- wave- like feature of the excitations, ft is noticed 
that the excitation of the spin-f part and that of the 
spin- 1/2 part are in phase for J' = — 4 but out of phase 
for J'=4. There is also clear difference between the two 
cases near the chain end. 
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FIG. 3. The excitation scheme (Sf) n - (Sf )i (for J' = -4) 
(Sf) n (for J'=4) on spin-1/2 chain (filled circle and square) 
and the spin-1 chain (open circle and square) for N=60. The 
circles are for the even i sites, while the squares are for the 
odd i sites. 

In Fig. ^, we show the correlation between neighboring 
site spins as a function of J'. The trends of these results 
are easy to understand by considering the proper limiting 
cases. When J' — > — oo, (Sfs%) approaches 1/6, show- 
ing that the cell composed by the nearest-neighbor (1,^) 
spin pair behaves like a S=3/2 spin; when J' — > oo, it 
approaches —1/3, the (l,i) spin pair is coupled as a spin 
singlet. The correlation (Sfsf +1 ) decreases monotoni- 
cally to a constant when J' — ► — oo, whereas when J' is 
positive, it first increases and reaches its maxim at J' ~ 1 
and then it decreases to zero for very large J' — > oo. 
The two intrachain correlation functions (SfSf +1 ) and 
reduce their magnitude when J' move away from 



zero. For J' — > — oo, they will decrease to non-zero con- 
stants while for J' — > oo, they will decrease to zero. 




FIG. 4. The dependence of the four neighboring correlation 
functions on the interchain coupling J' for J 1 — 1.0. 

We believe there is a phase transition between the 
strongly ferromagnetic rung coupling region and the 
strongly antiferromagnetic rung coupling region. The ex- 
istence of the end free spin reveals the VBS topological 
order in the ferromagnetic rung coupling region. Such 
a topological order is absent in antiferromagnetic rung 
coupling regies. Therefore a quantum phase transition 
should occur B3 We will further explore this phase tran- 
sition in a future study. 



III. MAGNETIZATION PLATEAU 

For a gapless system, the conformal field theorjEB 
shows that the upper and lower critical magnetic field 
H + and H ~ where the ground state of the system under 
magnetic field has magnetization m = M/N iaJ, 

H+ = lim [E(N, M + 1) - E(N, M)] 



N 

H- - lim [E(N, M) - E(N, M - 1)]. 

N-^-oo 



(4) 
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Here E(N, I) is the lowest energy in the subspace of an 
N-cell system with the z component of the total spin 
Ei(#i + sf) = /. If H+ - H~ - H m , there is no 
plateau at magnetization m in the thermodynamic limit 
and when the applied magnetic field is H m , the system 
has magnetization m. If and H ~ are not equal, there 
will be a magnetization plateau at m with a width D m = 
H^ —H^. By calculating the width of a length- N system: 

D m (N) = E(N, M + 1) + E(N, M — 1) — 2E(N, M), 

(5) 

D m , the plateau width at magnetization m can be ob- 
tained by considering the infinite length limit 
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D m = lim D m {N). 

N~* oo 



(6) 



According to Eq. (Q), Hamiltonian (||) may have a 
plateau at magnetization per cell m=l/2. Using DMRG, 
we have calculated Di(N) of this system and obtained 
two different plateau phases. 



A. J'=0 limit 

Before presenting the DMRG results, we examine the 
limiting case of J' = 0. In this case, the mixed spin 
ladder described by Hamiltonian (|^) decouples into in- 
dividual spin-1/2 and spin-1 chains. When a magnetic 
field is applied to the decoupled system, the spin-1 chain 
will not be magnetized until the field is larger than the 
critical field iff = J2As = i, where As=i is the Haldane 
gap of the pure spin-1 chain; while the gapless spin- 
1/2 chain will begin magnetization as soon as the mag- 
netic field is nonzero and reach its saturate magnetiza- 
tion at Hi = 2Ji. If Hf < HI, when the spin-1/2 

2 2 

chain has been saturated, the spin-1 chain has not been 
magnetized. Consequently, there will be a magnetization 
plateau between Hi and H\ , with a magnetization per 

2 

cell (S,s) of m=l/2. This plateau will happen when Ji 
is small and J 2 is large. Since we fix Ji + J 2 = 2, we 
can obtain a critical point for J\ where this plateau is 
diminished by taking 



^2^5=1 — Hi — iJ-j 3 = 2J\. 



(7) 



We have = 2+As \ • Haldane gap of the spin-1 
chain A# = i has been calculated with high accuracvEilE3, 
A s= i - 0.41, so Jf re 0.34. When J x < Jf, there is a 
plateau at m=l/2 with a width of D 1 / 2 = 2Ag = i — (2 + 
A S=1 )J X . 



B. DMRG results 

The energy levels E(N,N/2) E(N,N/2 + 1) and 
E(N,N/2 — 1) of Hamiltonian (||) are calculated using 
DMRG for different parameters J\ and J'. Because of 
the open boundary condition we use in.palculation, for 
Ji < Jj and small J' the edge statesEZlES will also ex- 
ist even at magnetization m=l/2. As a result, the low- 
est state in the N/2 + 1 sector may be the edge state. 
Considering this, we target at least two states for the 
E{N,N/2 + 1) sectors, and the plateau behavior corre- 
sponds to the second state of this sector E 2 (N, N/2 + 1). 

We check the accuracy of our results by examining the 
case of J 1 — and J' = 0, where the largest numerical 
error is expected for the DMRG calculation. In this case, 
the decoupled spin-1/2 sites are free spins. The m=l/2 
plateau begins at zero and lasts to the Haldane gap of 
the spin-1 chain, J2A s= i ~ 0.82. Our DMRG calculation 
gives E 2 {N,N/2 + 1) - E(N,N/2) ~ 0.832. 



When J' deviates from zero while J\ still is zero, the 
width of the plateau will decrease. We show 



Ht = E(N, N/2 -I 

2 a 

H+=E 2 (N,N/2 

2 

HI = E(N,N/2) 



l)-E(N,N/2), 
\-l)-E{N,N/2), 
- E(N, N/2 - 1) 



(8) 



at Ji = and J' = 0.6 in Figjol (a). Clearly Hf and i3T 

— 2 a 2 

converge to the same value in the thermodynamic limit. 

The plateau width is the difference otHf, and H^ , which 

2 " 2 

is about 0.42, smaller than that of the J'=0 case. For 
J' <0, the width decreases faster than in the J' >0 case. 
At critical points J' c _(Ji = 0) ~ 0.85 and J' c+ (Ji — 
0) ~ 2.2, the width decreases to zero and the plateau 
vanishes. In these cases, the width of the plateau can 
also be obtained from the difference of the E(N, N/2 + 1) 
and E 2 (N,N/2 + 1), D m = lim A r^ 00 [S 2 (A^, A^/2 + 1) - 
E(N,N/2 + 1)]. This definition gives D m with higher 
accuracy. The dependence of the width of the plateau on 
J' for Ji = is shown in Fig|| (b). 
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FIG. 5. Energy difference ((a) and (c)) Hf (filled cir- 

2 a 

cle), Hi (open circle), iff (filled square) at m—1/2 and the 

2 2° 

plateau width ((b) and (d)) in the two plateau phases. 

For larger J' > J' c , we find another plateau phase. The 
values of Hf , fft , and H~[ at J\ = and J' = 4.0 are 

shown in Fig.| (c). For this larger J', E 2 (N,N/2 + 1) is 
not an edge state. The non-zero magnetization plateau 
is obvious. The width of the magnetization plateau for 
Ji=0 and J' > J'c^h = 0) ~ 2.8 is shown in Fig. | (d). 

We have carried out systematic DMRG calculations 
to determine the magnetization plateau phase diagram 
shown in Fig. ^|. An examination of the calculated re- 
sults indicate that there arc two types of plateau phases 
corresponding to different magnetic structures in the lad- 
der. In phase I, the magnetization is due to the partially 
polarized s=l/2 spins. In phase II, the spin-1/2 sites 
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and spin-1 sites are strongly coupled. In phase III, where 
no plateau is found, the system is in a phase similar to 
that of a pure spin-3/2 chain. The magnetic features of 
the two plateau phases in Fig. || can be clearly seen in 
Fig. which shows the magnetization distribution on 
the spin- 1/2 chain defined as 



mi = — y^(s- 

2 ft ^— ' 



(9) 



The magnetization distribution on the spin-1 chain can 
be derived from mi = m — mi, where m is the total 

2 

magnetization of the ladder. 
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FIG. 6. The magnetization plateau phase diagram, the 
open circles are the critical points obtained by DMRG calcu- 
lations. 
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FIG. 7. Spin-1/2 chain magnetization mi at m=l/2 for 
Ji = (filled circle), 0.5 (open circle), 1.0 (filled square), 1.5 
(open square), and 2.0 (filled diamond). The solid lines are 
the spin-1/2 chain magnetization at J' — > ±oo limit. 

To understand the trends shown in Fig. |^, it is 
instructive to examine the two limiting cases. When 
J' — > — oo, mi — * 1/6 and mi = m — mi — ► 1/3, 
and the magnetizations on the two chains are half of 
their saturate values. In this limit, each cell (1,1/2) is 
in |3/2,l/2) = ^(V2|l,l;l/2,-l/2) + |l,0;l/2,l/2)) 

state when m=l/2, and the 1/2-spin and the cor- 
responding 1-spin bond together to form a 3/2-spin. 



When J' — ► oo, mi — > —1/6, that is, the orienta- 
tion of the spin-1/2 partition will be antiparallel to 
the applied magnetic field. Here each (1,1/2) cell is 
in 11/2,1/2) = 1/2,-1/2) - V2|l, 0; 1/2, 1/2)) 

state when m=l/2. 

We also have studied the magnetization process by ex- 
amining short chains. We show the typical plateau and 
no-plateau process in Fig.|](a) and Figjl(b) with Ji=1.0, 
J' = —1.0 and 3.0 respectively. For — 3.0 (phase II 
in Fig. |6|), the existence of the magnetization plateau at 
m = 1 is clear. When the chain length grows, hy 2 be- 
comes larger and /i7y 2 becomes smaller, but they converge 
to different values for the infinite chain. For J' = —1.0 
(phase III in Fig. ^), both hf^ 2 and h^M increase with 

chain length, but h~±, 2 increase faster. They eventually 
converge to the same value hin, yielding no plateau in 
this case. 




FIG. 8. Magnetization process of short chains with N=8, 
9, 10, 11, and 12 for Ji=1.0 and (a) J' = -1.0, (b) J'=3.0. 
The magnetization plateau at m=l/2 is obvious for J'=3.0. 

For J\ = J% = 1 with magnetization 1/2, it is known 
by perturbation that there is no spin gap for strong fer- 
romagnetic rung coupling region, and there is a finite 
spin gap for strong antiferromagnetic rung coupling re- 
gion. We have reported calculations for the intermediate 
rung coupling region by DMRG method in this section. 
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We find that there is only one quantum phase transi- 
tion at which the spin gap appears when we increase 
the rung coupling J' from a negative magnitude to a big 
and positive one. Beside this quantum phase transition 
between plateau phase II and plateauless phase III, the 
phase transition between plateau phase I and plateauless 
phase III is also a continuum phase transition due to spin 
gap closing. A more precise understanding of the gap 
generation and magnetization plateaus in the phase dia- 
gram found in this section can be given by bosonization 
discussions. We expect that further bosonization stud- 
ies will provide an interpretation of the plateau phase II 
beyond the obvious origin of its plateau. 

IV. CONCLUSION AND DISCUSSION 

We have calculated the low-energy properties and 
magnetization plateau behavior of a 2-leg mixed spin 
(S',s)=(l,l/2) ladder. From the low-energy properties, 
we conclude that for strong ferromagnetic interchain cou- 
pling, the 2-leg mixed spin ladder is in the same phase 
as the spin-3/2 chain. There is a magnetization plateau 
for large enough positive J' but no plateau is found for 
large negative J'. This suggests that at large J' > 0, the 
ladder would behave differently from the spin-3/2 chain. 

We have obtained the magnetization plateau phase di- 
agram of Hamiltonian (||) and identified two types of 
plateau phases. In the partially polarized spin-1/2 phase, 
the total magnetization comes mainly from the spin- 
1/2 chain, while the spin-1 chain is only polarized very 
weakly. In another plateau phase, because of the anti- 
ferromagnetic interchain coupling, the orientation of the 
spin-1/2 chain is antiparallel to the magnetic field, and 
the spin-1/2 chain and the spin-1 chain are coupled to- 
gether. 

Finally, we discuss some qualitative trends for gen- 
eral mixed spin (S,s) ladders. The basic structure of 
the phase diagram shown in Fig. || is expected to re- 
main valid. For higher spins (S,s), the non-plateau region 
phase III for (1,1/2) may also become a plateau phase. If 
S and s are both half-integers, there will be no partially 
polarized plateau phase. If S and s are one integer and 
one half-integer, then the low energy spectrum is gap- 
less and there is no plateau at m=0. In other cases, the 
low-energy spectrum is gaped, with a m=0 plateau. 
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